Steady streaming between two vibrating planes 
at high Reynolds numbers 

Konstantin Ilin and Andrey Morgulis 

Abstract 

We consider incompressible flows between two transversely vibrating solid walls and construct 
an asymptotic expansion of solutions of the Navier-Stokes equations in the limit when both the 
amplitude of vibrations and the thickness of the Stokes layer are small and have the same order of 
magnitude. Our asymptotic expansion is valid up to the flow boundary. In particular, we derive 
equations and boundary conditions, for the averaged flow. In the leading order, the averaged flow 
is described by the stationary Navier-Stokes equations with an additional term which contains 
the leading-order Stokes drift velocity. In a slightly different context (for a flow induced by an 
oscillating conservative body force), the same equations had been derived earlier by Riley [3]. The 
general theory is applied to two particular examples of steady streaming induced by transverse 
vibrations of the walls in the form of standing and travelling plane waves. In particular, in the 
case of waves travelling in the same direction, the induced flow is plane-parallel and the Lagrangian 
velocity profile can be computed analytically. This example may be viewed as an extension of the 
theory of peristaltic pumping to the case of high Reynolds numbers. 
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1 Introduction 

O , 

^ \ In this paper we study oscillating flows of a viscous incompressible fluid between two solid walls 
produced by the transverse vibrations of the walls. It is well-known that high-frequency oscillations of 
the boundary of a domain occupied by a viscous fluid generate not only an oscillating flow but also a 
(relatively) weak steady flow, which is usually called the steady streaming (see, e.g., the review papers 
[I] and [2, 3j). Recently such flows attracted considerable attention in the context of application of 
■ steady streaming to micro-mixing [H [6] and to drag reduction in channel flows [7j . 
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The basic parameters of the problem are the inverse Strouhal number e and the streaming Reynolds 



number R s , defined as 

V* V* 2 

* = \ Rs = \ W 
q ■ cod cov* 

where Vq is the amplitude of the velocity of the oscillating walls, d is the mean distance between 
the walls, u is the angular frequency of the vibrations and v* is the kinematic viscosity of the fluid. 
Parameter e measures the ratio of the amplitude of the displacement of the vibrating wall to the mean 
dist ance between the walls and is assumed to be small: e < 1. Note that the relation between the 
standard Reynolds number R based on V^* and d and the streaming Reynolds number R s is given by 
R = V^d/v* = R s /e, so that R s ~ 1 corresponds to R S> 1 when e is small. The present paper deals 
with the flow regimes with R s ~ 1 . This implies that the amplitude of the wall displacement is of the 
same order of magnitude as the thickness of the Stokes layer. 

Steady streaming at R s ~ 1 (and even for large R s ) induced by translational oscillations of a rigid 
body in a viscous fluid had been studied by many authors (see [El EJ [TUJ, EH E2J H3J E]). In the case 
of a body in an infinite fluid, the problem can be reduced to an equivalent problem about a fixed 
body placed in an oscillating flow, which considerably simplifies the analysis. A related problem of 
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steady streaming and mass transport produced by water waves had been treated in [171118]) where flow 
regimes with large R s were considered. Previous studies of flows produced by transverse oscillations of 
solid walls had been mostly focused on the problem of peristaltic pumping in channels and pipes under 
the assumption of low Reynolds numbers (R -C 1) and small amplitude-to-wavelength ratio (see, e.g., 
[Pol \W\ [6]). Although there are a few paper where the case of large Reynolds numbers (R 3> 1) had 
been considered [UG2E], these studies still correspond to R s <C 1. As far as we know, the flow regimes 
with R s ~ 1 have not been treated before. 

The aim of this paper is to construct an asymptotic expansion of the solution of the Navier-Stokes 
equation in the limit e — > and R s ~ 1. The procedure includes the derivation of the equations and 
boundary conditions for a steady component of the flow (steady streaming) that persists everywhere 
except thin layers near the vibrating walls (whose thickness is of the same order as the amplitude of 
vibrations). Since the thickness of the Stokes layer 5 = W u* Juj and the amplitude of the displacements 
of the walls a are of the same order of magnitude, the boundary conditions on the moving walls cannot 
be transferred to the fixed mean positions of the walls. This is what makes the problem difficult. 

To obtain an asymptotic expansion, we employ the Vishik-Lyusternik method (see, e.g., |19t l2L)j) rather 
than the standard method of matched asymptotic expansions. In comparison with the latter, the 
Vishik-Lyusternik method does not require the procedure of matching the inner and outer expansions 
and the boundary layer part of the expansion satisfies the condition of decay at infinity (in boundary 
layer variable) in all orders of the expansion (this is not so in the method of matched asymptotic 
expansions where the boundary layer part usually does not decay and may even grow at infinity). The 
Vishik-Lyusternik method had been used to study viscous boundary layers at a fixed impermeable 
boundary by Chudov [21]. Recently, it has been applied to viscous boundary layers in high Reynolds 
number flows through a fixed domain with an inlet and an outlet [22], to viscous flows in a half-plane 
produced by tangential vibrations on its boundary [23J and to the steady streaming between two 
cylinders [23]. A similar technique had been used in |25j to construct an asymptotic expansion in a 
problem of vibrational convection. 

The outline of the paper is as follows. In Section 2, we formulate the mathematical problem. In Section 
3, the asymptotic equations and boundary conditions are derived. Section 4 outlines the construction 
of the asymptotic solution. In Section 5, we consider examples of the steady streaming induced by 
vibrations of the walls in the form of standing or travelling waves. Finally, discussion of results and 
conclusions are presented in Section 6. 



2 Formulation of the problem 



We consider a three-dimensional viscous incompressible flow between two parallel walls produced by 
normal, periodic (in time) vibrations of the walls. Let x* = (x* , y* , z* ) be the Cartesian coordinates, 
t* the time, v* = (u* ,v* ,w*) the velocity of the fluid, p* the pressure, p the fluid density, v* the 
kinematic viscosity and d the mean distance between the walls. We introduce the non-dimensional 
quantities 

x* v* p* 

T = Ult , X = — , V = , p 



a uj a padu z 

where co and a are the angular frequency and the amplitude of the vibrations. In these variables, the 
Navier-Stokes equations take the form 

v T + e(v-V)v = -Vp + f 2 i/Vv, V-v = 0. (2.1) 

where 

a v* 
e = v 



d coa 
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Parameter e is the ratio of the amplitude of vibrations to the mean distance between the walls, v = i?" 1 
where R s is the Reynolds number based on the amplitudes of the velocity and the normal displacement 
of the vibrating boundary (the streaming Reynolds number). The walls are described by the equations 

z = ef(x,y,r) and z = 1 + eg(x,y,r), (2.2) 

where f(x,y,r) and g(x,y,r) are given functions that are 2-7r-periodic in r and have zero mean value, 
i.e. 

f(x,y) = J f(x,y,T)dr = 0, g(x, y) = j g(x, y, r) dr = 0. (2.3) 

o o 

We also assume that f(x,y,r) and g(x,y,r) are periodic both in x and in y with periods L x and L y , 
respectively. Boundary conditions for the velocity at the walls are the standard no-slip conditions: 

v\ z=ef = f T (x,y,r)e z , v \ z=1+eg = g T (x, y, r) e z . (2.4) 

The incompressibility of the fluid implies that functions / and g must satisfy an additional condition 
which, in the periodic case, is given by 

f(x,y,r)dxdy= g(x,y,r)dxdy. (2.5) 



oo oo 

The non-dimensional equations and boundary conditions depend only on two parameters e and v. In 
what follows we are interested in the asymptotic behaviour of periodic solutions of Eqs. (12. ip . (|2.4h 
in the limit when e — > and v = O(l). In other words, we consider the situation where the amplitude 
of the vibration is small and of the same order as the thickness of the Stokes layer. 

It is convenient to separate vertical and horizontal components of the velocity as follows: 

v = u + we z , u = ue x + ve y , 

i.e. u is the projection of the velocity onto the horizontal plane (the xy plane). We will also use the 
notation: r = (x, y) and x = (x, y, z). 



We seek a solution of (|2.ip . (|2.4p in the form 

u = u r (r, z, t, e) + u a (r, f , r, e) + u 6 (r, rj, r, e), 
w = w r (r,z,T,e) + ew a (r,£ i ,T, e) + ew b (r,i],T, e), 

p = p r (r, z, r, e) + p a (r, £, r, e) + p b (r, 77, r, e). (2.6) 

Here £ = z/e and 77 = (1 — z)/e are the boundary layer variables; all functions are assumed to be 
periodic in r, x and y with periods 2ir, L\ and L2, respectively. Functions u r , vf , p r represent a regular 
expansion of the solution in power series in e (an outer solution), and {u a ,w a ,p a } and {u b ,w b ,p b } 
correspond to boundary layer corrections (inner solutions) to this regular expansion. We assume that 
the boundary layer parts of the expansion rapidly decay outside thin boundary layers, namely: 

u a ,w a ,p a = o(£~ s ) as £ ->■ 00 and u b , w b ,p b = o(r]~ s ) as 77 ->• 00 (2.7) 

for every s > 0. In other words, we require that at a distance of order unity from the wall, the boundary 
layer corrections are smaller than any power of e. This assumption will be verified a posteriori. We 
begin with the regular part of the expansion. 
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3 Asymptotic expansion 

3.1 Regular part of the expansion 

Let 

V =vl + ev\ + e\ r 2 + f =p r + ep r 1 + e 2 p r 2 + (3.1) 

where v r = u r + w r e z and v r = u£ + w r k e z (k = 0, 1, 2, . . . ). The successive approximations v£, p r k 
(k = 0, 1, 2, ... ) satisfy the equations: 

d T vi = -Vpl + Fi(x,T), V-vg = 0, (3.2) 

where F{J(x,t) = 0, F;(x,t) = -(vg • V)vg and 

fc-i 

Pj(x, r) = - £(v$ • V)v^_ i _ 1 + ,V 2 vL 2 (A; = 2,3...). 
In what follows, we will use the following notation: for any 2-7r-periodic function /(r), 

/(r) =/ + /», /"= ^ | /(r)dr (3.3) 

o 

where / is the mean value of /(r) and, by definition, f(r) = f(r) — f is the oscillating part of /(r) 
that has zero mean value. 

3.1.1 Leading-order equations 

Consider Eqs. (|3.2p for = 0. We seek a solution Vq which is periodic in r. It can be written as 

v£ = % + v , v5 = V0 O (3.4) 

where 0o has zero mean value and is the solution of the boundary value problem 

V 2 </>o = 0, 4> (x + L x ,y,z) = 4> (x,y,z), 4>q{x, y + L y , z) = 4>q(x, y, z), 

4>0z \ z=0 = f T (x,y,r), cf> 0z \ z=l = g T (x,y,r). (3.5) 

Boundary conditions for (pQ Z at z = and z = 1 will be justified later. On averaging Eqs. (|3.2p for 
= 1, we obtain 

(v5-V)v5 + (v5-V)v5 = -VpI, V-v5 = 0. (3.6) 



Since, according to (|3.4j) . Vq is irrotational, we have: (vg ■ V)vq = (V (fro ■ V)V 4>o = V(| V^o| 2 /2) . 
Therefore, we can rewrite (|3.6[) as 

(vg • V)v5 = -Vn , V • v5 = 0, (3.7) 



where Ho = p[ + |V<^o| 2 /2- Equations p.7p represent the time-independent Euler equations that de- 
scribe steady flows of an inviscid incompressible fluid. It will be shown later that boundary conditions 
for Vq are 

v5| 2=Q =0, v5| 2=1 =0. (3.8) 

Thus, Vq(x) is a steady solution of the Euler equations (|3.7|) . which is periodic in x and y and subject 
to zero boundary condition at z = and z = 1. We choose zero solution: 

vS = (3.9) 

(evidently, it satisfies the Euler equations as well as the boundary conditions). This choice implies 
that there is no steady streaming in the leading order of the expansion. 
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3.1.2 First-order equations 

Separating the oscillatory part of Eqs. f|3.2|) for k = 1, we find that 



SrvJ = -vp;-(v5-v)v5, v-vi = o, (3.10) 

Using (j3.9[) and the fact that Vq is irrotational, ()3.10|) can be written in the form 

d T v[ = -VQi, V-vJ = 0, (3.11) 

where Q t = p\ + (Vcfo) 2 / 2 - It follows from ([3"XTj) that 

vj = vj + vi, v^ = V0i (3.12) 

where 0i has zero mean value and is the solution of the problem 

V 2 </>i = 0, 4>i(x + L x ,y,z) = <j)x(x,y,z), <j>i(x, y + L y , z) = 0i(x,y, z), 

<t>iz \ z=0 = ai(x,y,r), <j) lx \ z=1 = h(x, y, r), (3.13) 

where functions a%(x, y, r) and bi(x,y,r) will be defined later. 

On averaging equation for V3 and using (13.4p and (13.9P , we find that 

(vj • v)vj = - vm; + i/v 2 v^ -w r 2 x V0 O , v-vj = o. (3.14) 



where Ilg = P3 + (V</>i) 2 /2 + vJJ • V^o and u>2 = V x v^. It is clear from Eq. (I3.14f) that in order to 
obtain a closed equation for v£, we need to find v^. 

3.1.3 Second-order equations 

Consider Eqs. (|3.2p for k = 2. Averaging yields 



= -vf 2 - (vg • v)v; - (vs • v)vf, v • = 0. 

In view of (|3.4p . (|3.9j) and (|3.12p . the first of these reduces to the equation V (p r 2 + V0o ■ V0i) = 
that can be integrated to obtain 

$2 = — V(/>o • V(/>i + const. (3.15) 
The oscillatory part of (|3.2p for k = 2 gives us the equations for v£: 

v r 2r = -VQ 2 -(v5.VK-(v^V)v5, V-v£ = 0, (3.16) 

where Q2 = V2 + ■ V0i, and we have used (13. 4p . (13. 9p and (I3.12p . Taking cur/ of the first equation 
and using (|3.4p . we obtain 

w5r = [^5,«I] (3-17) 

where a>^ = V x and where [a, b] = (b • V)a — (a • V)b = V x (a x b) for any divergence- free vector 
fields a(x) and b(x). Now let £(x, r) be such that 

C T = vg, C = 0. (3.18) 

Then, it follows from (|3.17p that 

u, r 2 = [£M]- (3-19) 
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3.1.4 A closed system of equations for V][ 
In view of f|3. 18j) and (|3.19p . we have 

- G)\ x V^o = C T x *2 = C x [C, Q\\ = -C x [C T ,u{] = C x [QlXr], (3-20) 
where we have used the fact that /'(t)<7(t) = — /(t)</(t) for any periodic functions / and g. It 
follows from (pOOD that -£ r 2 x Vcfo = 5 (Cr x [C^i] + C x K,Cr]) • With the help of the identity 
ax [b, c] + c x [a, b] + b x [c, a] = V (a • (b x c)) (which is valid for any divergence-free a, b and c), 
this can be simplified to 

"*2 X V0 O = \ ICTTCT X S>\ + V (Q\ ■ (Cr X 0) • 

Finally, substituting the last formula into Eq. (|3.14p . we obtain 

(vj • V)vJ = -Vn 3 + i/V 2 v£ +Vx4 (3.21) 
where II3 = II| — Q\ ■ (C T x £) and 

V = |lC77<T- (3.22) 

When V is zero, (|3.2ip coincides with the stationary Navier-Stokes equations. In the two-dimensional 
case, (|3.22p reduces to the equation derived earlier in [3]. It was also noticed in [3] that V represents 
the Stokes drift velocity of fluid particles. 

3.2 Boundary layer equations 

In this subsection we will obtain asymptotic equations that describe boundary layers near the vibrating 
walls. 

Boundary layer at the bottom wall. To derive boundary layer equations near the bottom wall, we 
ignore u b , w b and p b , because they are supposed to be small relative to any power of e everywhere 
except a thin boundary layer near y = 1, and assume that 

u = u5 + ug + €(u; + u?) + ..., w = w r + e(w r 1 + w%) + ..., p = p r Q +pg + e(jp\ +pf) + . . . (3.23) 

We insert ()3.23j) into Eq. (j2.ip and take into account that u r k , w k , p r k (k = 0,1,...) satisfy the 
equations ()3.2p . Then we make the change of variables z = e £, expand every function of e £ in Taylor's 
series at e = and collect terms of the equal powers in e. This produces the following sequence of 
equations: 

« + ^| z=0 ^ + V||^-^|u^ = F^, d t pl = G a k , V r ug + %< = (3.24) 

for k = 0, 1, . . . Here V11 = e x d x + e y d y . Functions F k and G k are defined in terms of lower order 
approximations. In particular, Fg = 0, Gq = 0, G\ = and 

F? = - « + e<) U "of - K • V[|) L =0 ug - (ug • V,|) uS \ z=Q - (ug • V,|) ug - ufiufc. (3.25) 

Boundary layer at the upper wall A similar procedure leads to the following equations of the boundary 
layer near the upper wall: 

d T u b k -w r \ z=1 d v u b k + X7 llP b k -vd 2 v u b k = Fl d vP b k = Gl V v u b k -d v w b k = (3.26) 
for k = 0, 1 ... , where F b = 0, G b = 0, G\= and 

PS = K - V< z ) U < ~ K ■ V,|) L =1 u b - (ug ■ V„H \ z=l -(u b • V,|K + «■ (3-27) 
In accordance with ()2.7|) . we require that 

u t> w k,Pk = °(C S ) as C^oo and u|, w\, p\ = o(r/~ s ) as 77 -> 00 (3.28) 
for every s > and for each A; = 0, 1, ... . 
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3.3 Boundary conditions 

Now we substitute (|3.23[) in (|2.4p . expand all functions corresponding to the outer flow in Taylor's 
series at z = and collect terms of equal powers in e. This leads to the following boundary conditions 
at the bottom wall: 

u5| z=0 +ug| f=/ =0, w r \ z=Q =f T ; (3.29) 

h 

u£ \ Z=Q +u a k \ £=f = d>l-n Lo (^ = 1,2,...), (3-30) 

n=l 

wl \ z=0 = -wU | ?=/ -Y, J -t d>i- n \ z=0 (A; = 1,2,...). (3.31) 

71 = 1 

A similar procedure yields the following boundary conditions at the upper wall: 

11^+11(51^=0, wl\ z=1 =g T] (3.32) 

71 

»tf U +*4 - E |r a >£-n U (* = i, 2, . . . ), (3.33) 

n=l 
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\r,=-g ^ n , 
n=l 



E^T a >fc-nL=i (fc = l,2,...)- (3-34) 



Note that the boundary conditions for <j)Q Z at z = and z = 1 in problem (|3.5j) follow directly from 
(13391) and (l332j) . 

4 Construction of the asymptotic solution 
4.1 Leading order equations 

Oscillatory outer flow. In the leading order, the oscillatory outer flow is irrotational, and the velocity 
potential (f>o is a solution the boundary value problem (|3.5p (which is unique up to a constant). 

Boundary layer at the bottom wall. In the leading order, the boundary layer is described by Eqs. (|3.24p 
with k = 0. The decay condition (I3.28P for and the second equation (I3.24p imply that p$ = 0. 
Hence, the equation for reduces to 

ug r + / T ug £ = i/ug £{ . (4.1) 
Here we have used boundary condition (|3.29p for Wq. Let 

s = e-/(r,r) and ug(r, £, r) = Ug(r, s(r, £, r), r). 
Then Eq. (|4.ip simplifies to the standard heat equation for Ug(r, s,t): 

U a or = vU a ss- (4-2) 
Boundary conditions for Ug are the decay condition at infinity and the condition 

USLo=- u 5U> ( 4 -3) 

which follows from (|3.29p . Averaging Eq. (|4.2p . we obtain the equation Ug ss = 0. The only solution 
of this equation that satisfies the decay condition at infinity is zero solution Ug = 0. [Note that this 



7 



doesn't mean that u[j = 0, because Uq is averaged for fixed s, and u[j - for fixed £. This fact, however, 
is inessential for what follows.] 

A periodic (in r) solution of Eq. (|4.2p satisfying the decay condition at infinity and boundary condition 
(|4.3p and having zero mean can be found by standard methods (see, e.g., [26J). Once Ug is found, the 
normal velocity u>q is determined from the continuity equation (the third equation (|3.24p ) that can be 
written as 

V r Ug-Ug s -V||/ + ^ a s = (4.4) 

where function Wq(t,s,t) is defined by the relation Wq(t,^,t) = Wq (r, s(r, £,t), t). Integration of 
Eq. ([HID in s yields 

oo 

WS(t,s,t) = Ug-V||/ + V|| • Jv a (r,s',T)ds', (4.5) 

s 

where the constant of integration was chosen so as to guarantee that Wq — > as s — > oo. 
Boundary condition (|3.31|) for k = 1 is 

w\ \ z=0 = -fw r 0z \ z=0 ~ | e=/ or, equivalents, w\ \ z=0 = - f w r 0z \ z=Q - W$ \ s=Q . (4.6) 

This means that function a%(x,y,T) in the boundary value problem (|3. 13|) is given by 

a 1 (x,r) = -f^ z \ z=Q -W a \ s=0 . (4.7) 



Ul = vU b 0qq , (4.8) 



Boundary layer at the upper wall. Exactly the same analysis as in the case of the bottom wall leads 
to the heat equation 

Tib , , T , 

'Oqqi 

where q and XJ b are defined as 

q = r) + g(r,r) and Uo(r, 77, r) = Uo(r,g(r, 77, r),r). 
Boundary conditions for Ug are the decay condition at infinity and the condition 

U oU=- u 5Li> (4-9) 

that follows from (|3.32p . Equation (|4.8j) has a consequence that Uq 99 = 0. This, together with the 
condition of decay at infinity (|3.28p . implies that Uq = 0. 

Again, an oscillatory solution of Eq. (I4.8P can be found by standard methods. The normal velocity 
Wq is determined from the continuity equation (the third equation (|3.26|) ): 

00 

W 6 (r, q, r)=V b - Vyff ~V\\-J Ug(r, q' , r) dq', (4.10) 

where function W$ is defined by w^r, rj, r) = W$ (r, q(r, r], r), r). and where the constant of integration 
was chosen so as to guarantee that Wq — > as q — >■ 00. 

Boundary condition ()3.34p for k = can be written as 

W l\ z=1 =-gw r 0z \ z=1 -W b \ q=0 . (4.11) 
This means that function bi(x,y,r) in the boundary value problem (|3.13j) is given by 

b x {x,r) = -g^ z \ z=l -W b Q \ . (4.12) 
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Averaged outer flow. Averaging boundary conditions (|3.29p and (|3.32[) . we find that 



US | z=0 = " ug | ?=/ = -Ug \ s=0 = 0, W r \ Z=Q = f T = 0, 



I =- U [;| = — Ug I „=0, Wn\ .= g-=o. 

\z=l U \rj=— g u 1 1/=0 ' u lz=l J/ 



This justifies our earlier assumption (|3.8p and our conclusion f|3.9j) that there is no steady streaming 
in the leading order of the expansion. 

4.2 First-order equations 

Oscillatory outer flow. The oscillatory part of the (first-order) outer flow is determined by the bound- 
ary value problem ()3.13j) for the Laplace equation. Since functions ai(r,r) and &i(r, r), which appear 
in the boundary conditions for <pi z , are now known, problem (|3.13p can be solved, thus giving us 
vl(x,r). 

Boundary layer at the bottom wall Consider now the first-order boundary layer equations (Eqs. (|3.24p 
for k = 1). Again, the condition of decay at infinity for p\ and the second equation (|3.24p imply that 
p\ = 0, so that the first equation (I3.24p simplifies to 

Ulr+/r4"Kr F l' ( 4 - 13 ) 

We are looking for a solution of (|4.13p which satisfies the decay condition at infinity (|3.28p and 
boundary condition (I3.30P for k = 1. To find it, we again employ variable s = £ — /(r,r) and rewrite 
(j3~T3l) and (^30]) in the form 

U? T - v\J a lss = H?, (4.14) 

u?L =0 =-^L=o-/^L=o' ( 4 - 15 ) 

where Uf(r,s,r) is such that Ui(r, s(r, £, t),t) = uf(r,^,r) and where 

h? = -K + ( s + /)OL =0 us s -K-v)|^ ug + K-v/)| 2=0 ug s 

- (Ug • V) u5 | z=Q - (Ug • V) Ug + (Ug • V/) Ug s - WS Ug s . (4.16) 
Averaging (I4.14j) . we find that 

^U? ss = -H?(r, S ). (4.17) 

We integrate this equation two times, choosing constants of integration so as to satisfy the condition 
of decay at infinity in variable s. As a result, we have 

oo oo 

UJ = -i/" j-R a x (r,s")ds" ds'. (4.18) 

s s' 

Note that (I4.18P is the unique solution of Eq. (|4.17p that decays at infinity. The oscillatory part of 
Uf can also be found from Eqs. (|4.14p and (|4.15p . but we will not do it here as we are only interested 
in the averaged part of the flow. 



Boundary layer at the upper wall. A similar analysis leads to the equations 

>lqq 



U 6 lT -^U? =H?, (4.19) 



U i L=n= _u i L=i -9< z \ z =v ( 4 - 2 0) 



| g =0 1 \z=l » 1 2= i 
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where U*(r, g,r) is defined by the relation U^(r, q(rj,r), r) = u^(r, r/, r) and where 

H? = (w r 1 -(q-g)w r 0z )\ z=1 V b 0q -(u r .V)\ z=1 \J b -(u r .Vg)\ z=1 V b 0q 

U b o • V) uS | 2=1 - (ug ■ V) Uj - (Ug • V 5 ) + W* USg. (4.21) 



Averaging (|4.19f) . we obtain 
Integration of (I4.22p yields 



v\J b lqq = -n\{r,q). (4.22) 

oo oo 

V\ = -\f I H\(r, q")dq"dq'. (4.23) 

9 q' 

Here again the constants of integration are chosen so as to satisfy the condition of decay at infinity. 
Averaged outer flow. Averaging the boundary conditions (|3.30p and (|3.33[) for k = 1, we obtain 



l 2 =o-- U ?Lo-/ u U 2 =o' ( 4 - 24 ) 



u?U-»"5.U- ( 4 - 25 ) 



| 2 =1 ~±l 9 =o 



Once Uf and U\ are found (from Eqs. (jlTTHD and (ET231) ). Eqs. (I4T241 and (ET25D give us boundary 
conditions for u[. 

Boundary conditions for w\ are obtained by averaging Eqs. (|4.6p and (|4. 1 lj) : 



^L=o=-/^L =0 -ug-v/j s=0 , 



01 L=i= " 2«4 L=i - u o • v 2 1 9=0 • 

Here we have used Eqs. (|4.5|) and (|4.1U|) . With the help of (|4.3|) and (|4.9p . these can be transformed 
to 



^Lo=V||-(/V|]0o) 2=0 , (4.26) 



l*=l 



V||-( 5 V||0 Q )| z=1 . (4.27) 



Equations (|4.24p ~ (|4.27p represent boundary conditions for the first-order averaged equations for the 
outer flow (|33Tj) . 

Remark on boundary conditions (|4.26p . (|4.27p . In general case the right sides of Eqs. (|4,26p and (|4.27l) 
are nonzero. Nevertheless, the net flux of the averaged Eulerian velocity through any rectangle of 
periods in x and y always vanishes due to periodicity of /, g and 0o an d the form of the right sides 
of Eqs. (|4.26p . (|4.27[) . It can be shown (see Appendix A) that the z-component of the Stokes drift 
velocity at z = and z = 1 is equal to minus the right sides of Eqs. (|4.26p and (|4.27j) . so that the 
averaged Lagrangian velocity at z = and z = 1 vanishes. This means that at each point of the 
boundary the averaged mass flux of the fluid is exactly zero. 

5 Examples 

5.1 Example 1: standing waves 

Let the vibrations of the walls be in the form of standing waves: 

f(T,T) = R(r)h(r), g(T,T) = S(r)h(r), (5.1) 
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where R(r) and S(r) are given functions, periodic in x and y with periods L x and L y , respectively, 
and satisfying condition (|2.5p . and where h(r) is a given 27r-periodic function. In this case, the 
solution of problem (|3.5p has the form o (r, z,r) = <£(r,z) /i'(r) where <£(r, z) is the solution of the 
time-independent problem 



V 2 $ 



0, $(x + -C a; ,y,2;) = §(x,y,z), $(x, y + L y , z) = $(x, y, z), 



L=o= 

It is not difficult to show that 



<I>- 



I 2 = 1 



5(r). 



(5.2) 



fv r 0z \ z=0 =0, 9< z \z=i=°> V |r (/V||0o)| 2=o = 
Therefore, boundary conditions (|4.24p - (|4.27p simplify to 



0, V,| • (gVn 



u, 



I 2=0 



-U? 



s=0' 



u. 



lz=l 



-U 



1 9=0' 



I z=0 



0. 



1 I 2=1 



1 2=1 



0. 



(5.3) 



Thus, vibrations of the walls in the form of standing waves produce zero boundary conditions for the 
normal velocity in the averaged outer flow. 

It is also easy to deduce from (|3.22p that the Stokes drift velocity V is zero, so that Eq. (|3.2ip reduces 
to 



K • v) v\ 



-vn 3 + vv 2 v r 1 . 



Thus, in this example, the first-order steady outer flow is described by the Navier-Stokes equations 
with boundary conditions (15. 3D . 



In particular, if we consider vibrations of the walls in the form of plane standing waves 

f(x, y, t) = cos(kx) cos r, g(x, y,r) = a cos(fcx) cos r, 

where k = 2ir/L x is the wave number, a = ±1 (a = 1 corresponds to bending waves and a = —1 - to 
contraction/expansion waves, see Fig. [1]), then the boundary conditions (|5.3p will take the form 



'1U=0- 



- - k Q 2 (k) sin(2kx) e x , v 



where 



Q{k) 



1 \z=l 



1 — a cosh k 
sinh k 



3 

--kQ 2 (k) sm(2kx)e x , 



(5.4) 



The results of numerical solution of the Navier-Stokes equations subject to the boundary conditions 
(I5.4p are shown in Fig. [2j Qualitatively, the stream line picture is the same for both bending and 
contraction/expansion waves. The magnitudes of the averaged flow in these two cases are different 
for L x > 1. Let us measure the magnitude of the velocity field by max|v^| where the maximum is 
computed for all (x,z) such that < x < L x , < z < 1. Table 1 shows maxlvjl for various L x . 
Evidently, the maximum velocity decreases, when wavelength L x increases, much faster for bending 
waves (a = 1) than for contraction/expansion waves (a = —1). 



L x 


0.5 


1 


5 


max 


a = 1 
a = —1 


2.3562 
2.3562 


1.1693 
1.1869 


0.073073 
0.75974 



Table 1: max for standing waves and v — 1. 
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(a) a = 1 



(b) a = -1 



Figure 1: (a) bending waves; (b) contraction/expansion waves. 



(a) a = 1 , l_ x = 0.5 (b) a = 1 , L x = 1 (c) a = 1 , l_ x = 5 




Figure 2: Streamlines of the averaged velocity for v = 1 (standing waves). 
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(a) a=1 



(b) a=-1 





Figure 3: The averaged Eulerian velocity as function of k. 
5.2 Example 2: travelling waves 

Consider now vibrations of the walls in the form of plane travelling waves: 

/(r, r) = cos(fcr — r), g(r, r) = a cos(kx — f3r), 



(5.5) 



where k = 2n/L x , a = ±1 and (5 = ±1. If (3 = 1, both waves travel in the same direction, if (3 = — 1, 
the waves travel in opposite directions. Since both / and g do not depend on y, we will look for a 
solution which does not depend on y. Problem (|3 . 5f) simplifies to 



(5.6) 



4>0xx + 0022=0, (j) (x + L x ,z) = <j) (x,z), 
4>0z \ z=0 = sin(kx - r), (j) 0z \ z=1 = asin(kx - r). 



Consider first f3 = 1, i.e. both waves travel in the positive direction of the x axis. Boundary conditions 
(jTOD - flPTD become 



where 



__ Q =A(k), u\\ z=1 = A(k), w[\ z=Q =0, w[\ z=1 - 
k 



0. 



(5.7) 



A(k) 



4sinh 2 (fe) 



[5 + cosh 2 (/c) — 6a cosh(fc)] 



The Stokes drift velocity is given by 
V(r 



k [cosh(/c) — a] , r , 
V (z)e x , V (z) = 1 . \ ' x J cosh [k(2z - 1)] 



sinh (k) 

Now u)\ = fi(x, z)e y where £l(x, z) = dzU^ — d x w\, and Eq. (|3.2ip can be written as 

(v£ • V)vJ = -Vn 3 + uV 2 w r x + V (z)n e z . 
Equation (|5.9p subject to boundary conditions (|5.7p has the constant solution 

= A{k)e x . 



(5.8) 



(5.9) 



(5.10) 



The constant horizontal velocity as function of k is shown in Fig. [3l For bending waves, it is negative 
for small k and positive for large k (A(k) ~ k for both k <C 1 and A; ^ 1). For contraction/expansion 
waves, it is always positive and growing both when k — > and when k — > oo {A(k) ~ for A; 1 
and A{k) ~ A; for A; > 1). 
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(a) 0=1 (b) a=-1 




0.5 1 1.5 2 4 6 8 10 12 

Figure 4: The averaged Lagrangian velocity profiles for (3 = 1 (waves travelling in the same direction). 

Although Eq. (|5.9p also has solutions with a nonzero pressure gradient VH3 = Coe x (cq = const), 
we do not consider such solutions here, because this would amount to a modification of our original 
problem, allowing the presence of a weak 0(e 3 ) pressure gradient. Note also that flows with stronger 
pressure gradients would require a separate treatment. 

The averaged Lagrangian velocity is given by 

vf = v[ + V = V L (z)e x , V L (z) = A(k) + V (z). 

Its profiles for various values of k are shown in Fig. HI In the short wave limit, V L (z) is constant 
almost everywhere except narrow layers near the boundary where the velocity is much higher (see Fig. 
[5]). This can also be seen from the asymptotic formula 

V L = k(^ + e~ 2kz + e - 2k{1 "A + 0(ke~ k ), k -> 00. 



Now let /3 = —1. In this case, boundary conditions (j4.24j) f|4.27[) take the form 

u\ | 2=0 = B + (k,x), u\ \ z=1 = B~(k,x), w\ \ z=0 = C(k,x), w\ \ z=1 = C(k, x), (5.11) 



where 



B ± (k,x) 



cosh(/c) 

±1 + 6a s sin(2K2;j 

sinh (A;) 



, C(k, x) = a 



k sin(2A;2;) 
sinh(fc) 



(5.12) 



Thus in the case of waves travelling in opposite directions, we have nonzero boundary conditions for 
the normal velocity in the averaged outer flow. 



The Stokes drift velocity is given by 

V(r) = V 1 (z)e x + V 3 (x)e z , V x (z) 



smh[k(2z - 1)] 
~ k r T77^ ' V 3 (z) = -C(k,x), 



(5.13) 



where C(k,x) is given by (|5.12p . Note that Eq. (|5. 13|) implies that the normal component of the 
averaged Lagrangian velocity at the walls is zero: w = Wi + V3 = at z = and z = 1, which is in 
agreement with the general property discussed in Appendix A. 

Further, we have = Q(x,z)e y where, as before, Q(x, z) = d z u\ — d x w\, and Eq. (|3.2ip can be 
written as 

(9; • V)vJ = -Vn 3 + vV 2 v\ + n [V^z) e z - V 3 (x) e x ] . (5.14) 
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(a) L x = 1 (b) L x = 2 (c) L x = 3 




Figure 6: Streamlines of the averaged Lagrangian velocity for a = 1 and v = 1 (waves travelling in 
opposite directions). 
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Thus, the averaged Eulerian velocity can be found by solving Eq. (|5.14p subject to boundary conditions 
(|5,lip . Then the averaged Lagrangian velocity is obtained by adding the Stokes drift velocity given 
byEq. ([533]) . 

Equation (|5,14p supplemented with the incompressibility condition and boundary conditions (|5,lip 
were solved numerically. The streamlines of the averaged Lagrangian velocity v L for a = 1 and v = 1 
are shown in Fig. [6l The corresponding values of max |v£| are given in Table 2. In the case a = — 1, 
the same flow pattern is shifted by a quarter of the wavelength in the x direction. 





1 


2 


3 


5 


8 


15 


max 


0.80301 


0.60373 


0.71208 


1.0467 


1.6146 


2.9929 



Table 2 : max |vjj for (3 — —1, a — 1 and v — 1 (waves travelling in opposite directions). 



5.3 General periodic vibrations 

Here we will briefly describe how to treat the case of general periodic (both in time and space) 
transverse vibrations of the wall. We focus our attention on the travelling waves for which we will 
present an explicit solution. 

Let <r be the projection of £ (defined by Eq. (3.18)) onto the xy plane, i.e. 

c- = C-(C-e 2 )e 2 , (5.15) 

and let 

<r a = <r\ z=0 , ° b = °\ z=l - (5-16) 
Let us express boundary conditions (|4.24p . (|4,25p in terms of /, g, <r a and cr b . 

First, we note that since cr a and <r b are periodic functions of r with period 2tt, they can be written as 
Fourier series 

v = J2"n(r,z)e inT , <7° = 5>£(r)e^, <x fe = (r)e m \ (5.17) 

n n n 

Here and in what follows (unless explicitly specified otherwise) summation is performed over all inte- 
gers. Also, in Eq. (|5.17p . 

a (r,z) = 0, <xg(r) = 0, <x£(r) = 0. 



Now boundary conditions (|4.3p and (|4.9p can be written as 

UoUo^E™^' US| =-i£n*ne inT . (5.18) 



71 



Periodic solutions of Eqs. (14. 2j) and (|4.8p subject to boundary conditions (I5.18P and the condition of 
decay at infinity (in s and q respectively) are given by 

US = na a n e inT ~^ s , XJ b = -ij^ n& b n e inT -^ q , (5.19) 

71 71 

where 




Mn = \/7r 1 (l + isgn(n)). (5.20) 
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It is shown in Appendix B that boundary conditions for (given by (|4.24p and (|4.25p ) can be written 

as 



I 2=0 



-|«t«(V|| • a*) - \V\\ y£ \n\\a a n A - M Re [<*(V|| • <)] , (5.21) 

\ n / n 

u[\ z=1 =g7^\ z=1 -^(V r ^)-iV|| (j2\ n W & n\ 2 ) -EH Re [^ (V " - &b A ■ (5 ' 22) 

\ n / n 

From now on we will restrict our analysis to the two-dimensional case, where both / and g do not 
depend on y and = u\e x , a = ae x . Equations (|5.2ip and (15. 22ft simplify to 



l r i\ z=0 = fr<r z \ z=0 ~a-a--^d x (l>|K| 2 J > ( 5 - 23 ) 
1{\ Z=1 =W^\ Z=1 -\v^i-\d x (^\n\\a b n A . (5.24) 



For simplicity, we consider only vibrations of the walls in the form of plane waves of arbitrary shape 
travelling in the same direction^, namely, we assume that 

f(x,T) = F(kx-r), g(x,r) = aF(kx-r), (5.25) 

where F is an arbitrary 27r-periodic function of a single variable h = kx — r; k and a are the same 
as those defined in Section 5.2. Let F n be the Fourier coefficients of F. Then the solution of problem 

<h(z,h) =Y J ™b n {z)F n J nh (5.26) 



where 



Further, we have 



and 



w \ cosh[|fen|(l — z)\ — a cosh[|/cn|z] 
\kn\ sinh \kn\ 

a(z, h) = -ik ^ nb n (z)F n e inh (5.27) 



a a (h) = iak^2nQ n F n e inh , a b (h) = -tk^nQ^e 1 ^, Q„ = L-^g^M , ( 5 .28) 

It follows from (l5T28|) that 

ot n (x) = iknaQ n F n e inkx , a b _ n {x) = -iknQ n F n e inkx . (5.29) 
Substitution of (|5.27p - (|5.29p into (|5.23p . (|5.24p and tedious but elementary calculations yield (cf. Eq. 



=1 =A(fc), 



, x — 9, a ,9 5 + cosh 2 \kn\ — 6a cosh \kn\ 

^*) = *E»W 2 'J, a|M 

_ 2 1 f 1 2 5 + cosh2 (kn) -6a cosh(kn) 

- k l^ n 1^1 S inh 2 (fcn) • (5 ' 30j 

n=l v ' 



More general vibrations can be treated similarly. 
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Boundary conditions for w\ (given by (|4.26p . (|4.27p ) reduce to 

m =0 =0, ^| 2=1 =0. (5.31) 
Similar calculations of the Stokes drift velocity yield 

V = V (z)e x 



where (cf. Eq. (15.8 



V = 2k y n 2 \F n \ 2 COSh |fc 2 re| ~ - cosh [\kn\(2z - 1 
^-^ sinh \kn\ 



n 
oo 



= 4A;V n 2 |Fn| 2/ : cosh [fcre(2z - 1)] . (5.32) 

z — ' sinh (kn) 

n=i y ' 

Boundary conditions (15.30P and (|5.3ip results in a constant solution of Eq. (15. 9ft : = A(k)e x , which 
is a generalisation of the solution (|5.10p obtained in Section 5.2. The averaged Lagrangian velocity is 

V L = V L e x , V L = A + V (z). 

Let us compute the total volume flux of the fluid through the channel (per unit length in y direction) . 
We have 

Q= V L dz = A+ I V (z)dz = yV|F n | 2 R n {k) (5.33) 
Jo Jo n=1 

where 

R n = ^ 5 + cosh 2 (A;n) - 6a cosh(/cn) + 4 [cosh(fen) - a] ^ — - . (5.34) 

sinh z (/cn) \ kn J 

Optimal shape of the wave. Now we can answer the question on what shape of the travelling waves 
leads to the the most efficient generation of the averaged flow. Mathematically, this can formulated as 
the following optimization problem: maximize the total volume flux Q on the set of all F n such that 

oo 

^V|F„| 2 = 1. (5.35) 

n=l 

Let m € N be such that 

R m = maxR n . (5.36) 

nGN 

It is not difficult to show that Q attains its maximum for 

e iX 

F n = -ri5\rim (n = ±1, ±2, . . . ), 
\n\ 

where A is any real number. For contraction/expansion waves (a = —1), the sequence {R n } is 
monotonically decreasing for all k > 0, so that the maximum is attained at n = 1. Thus, the most 
efficient deformation of the walls is given by (cf. (15.6|) ) 

f(x,r) = 2cos(fcr -r + A), g{x,r) = -/(x,r). 

In the case of bending waves, the mode for which the maximum total flux is attained depends on k 
as shown in Fig. [71 One can see that m is equal to 1 for k > 2.9 (i.e. for sufficiently short waves) and 
increases when k decreases. The most efficient deformation has the form: 

f(x, r) = g(x, r) = 2 cos[m(kx — r) + A]/m. 
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Figure 7: to, defined by (|5.36p . as a function of k for bending waves (a = 1) 



6 Discussion 

We have considered incompressible flows between two transversely vibrating solid walls and proposed a 
general procedure for constructing an asymptotic expansion of solutions of the Navier-Stokes equations 
in the limit when both the amplitude of vibrations and the thickness of the Stokes layer are small and 
have the same order of magnitude. The procedure is based on the Vishik-Lyusternik method and, 
in principle, can be used to construct as many terms of the expansion as necessary. In the leading 
order, the averaged flow is described by the stationary Navier-Stokes equations with an additional 
term which contains the leading-order Stokes drift velocity. In a slightly different context (for a flow 
induced by an oscillating conservative body force), the same equations had been derived earlier by 
Riley [3]. 

The general theory has been applied to two particular examples of steady streaming induced by 
vibrations of the walls in the form of standing and travelling plane waves. In the case of plane 
standing waves (bending or contraction/expansion waves), the averaged flow has the form of a double 
array of vortices (see Fig. [2]). For short waves the vortices are concentrated near the walls, while long 
waves produce vortices that fill the entire channel. The latter pattern is similar to flows produced 
by long standing waves at small Reynolds numbers [6]. For both bending and contraction/expansion 
waves, the flow patterns are qualitatively the same. However, the intensity of the averaged flow is 
different. Although it is of similar magnitude for both waves when the wavelength is small, for long 
waves, the maximum velocity decreases much faster for bending waves than for contraction/expansion 
waves, which is in agreement with the classical theory of peristaltic pumping at low Reynolds numbers 
(see, e.g., [15]). 

If vibrations of the wall have the form of plane harmonic waves which travel in the same direction 
(this leads to a bending wave if the phase shift between these travelling waves is zero, and to a 
contraction/expansion wave if the phase shift is w), the induced steady flow is a unidirectional two- 
dimensional flow. The fortunate thing is that, in this case, the asymptotic equations can be solved 
analytically. The Lagrangian velocity profile is symmetric relative to the channel central axis, with 
a minimum velocity at the centre and maxima at the walls. In the short wave limit, the Lagrangian 
velocity profiles for both bending and contraction/expansion waves are the same: the velocity is nearly 
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constant across the channel except for narrow layers near the walls where velocity is much higher. For 
moderate and long waves, the velocity variation across the channel is much smaller, and the magnitude 
of the velocity rapidly grows for the contraction/expansion wave and decays for the bending wave when 
the wavelength increases. This example may be viewed as an extension of the theory of peristaltic 
pumping to the case of high Reynolds numbers. 

If vibrations of the walls have the form of plane harmonic waves traveling in opposite directions, the 
averaged flow depends on the wavelength in a more complicated way. For short waves the averaged 
flow is a superposition of a shear flow with a linear velocity profile and a periodic array of weak vortices 
('cat's eyes') in the center of the channel. When the wavelength increases, the vortices first grow in 
size and magnitude. Then each of them splits into a pair of vortices of the same sign, and at the same 
time new vortices appear near the walls and move towards the centre of the channel. Eventually, for 
long waves, there are two arrays of alternating vortices (similar to what we had in the case of standing 
waves) . 

As an example of general periodic vibrations, we have considered travelling contraction/expansion and 
bending waves of an arbitrary shape. Similarly to the case of a harmonic vibration, these produce 
unidirectional two-dimensional flows. The averaged Lagrangian velocity is written as an infinite series 
in term of Fourier coefficients of the deformations of the walls. The natural question that appears here 
is: what deformations of the walls lead to a maximum total volume flux through the channel. It turns 
out that the optimal shape always corresponds to a single harmonic. For a contraction/expansion 
wave, it is always the first harmonic (with m = 1). For a bending wave, the optimal harmonic is 
uniquely determined by the parameter k = 2ir/L x . For example, if k = 1, then m = 4. It should 
be noted here that a contraction/expansion wave always produces a higher total volume flux then a 
bending wave. 

There are many open problems in this area. First, it is not quite clear how the present theory can be 
extended to the case of high R s . Second, we did not make any assumption about the characteristic 
length scale of vibrations in the horizontal direction. Nevertheless, the examples show that the theory 
does not work when the horizontal length scale is much smaller than the mean distance between the 
walls. In particular, Figures 2 and 5 shows that for L x < 1 the flows in Examples 1 and 2 become 
concentrated near the walls. This suggests that in the limit of short waves, one can try to construct a 
theory with double boundary layers (similar to what had been done by Stuart [8] and Riley [9] for an 
oscillating cylinder at R s 3> 1). Third, it would be interesting to investigate the effect of an external 
mean flow (e.g., a flow produced by a mean pressure gradient) on the steady steaming between two 
walls and vice versa. This problem has an additional parameter - the ratio of the characteristic velocity 
of the mean flow to the amplitude of the velocity of the vibrating walls and, therefore, one can expect 
a variety of different flow regimes. All these are problems for a future investigation. 

Acknowledgments. This work was initiated during a short visit of Andrey Morgulis to the University 
of York under the University of York Research Development Visit scheme. 
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7 Appendix A 

Here we will show that 



^3 U= "V|| • (/V,|0o) | 2=0 , V 3 U= -V,| • (gV^o) \ z=v (A.l) 
where V3 is the z-component of the Stokes drift velocity V given by Eq. (|3.22[) . 
We have 

1 



^3 = V-e, = -e 2 .[C T ,C] = -(C-V)(C-e z ) = -V-[C(C-e z )] 



= -d z [(C T -e z )(C-e z )]-V v [u r (C-e z )] = -V v [u r (C-e z )]. (A.2) 

It follows from (|3.18p and (|3.5p that £ • e z \ z= Q = f and £ • e 2 | 2= i = g. These and (|A.2p imply the 
relations 

^L = o=-V| r CWUo and V 3 \ z=1 =-V v (g^)\ z=0 
that are equivalent to (jA.ip . 

8 Appendix B 

Here we will briefly describe the derivation of Eq. (|5.2ip . Equation (|5.22p can be derived in exactly 
the same way. 

First, we use (14. 9p and (14. lip to rewrite Eq. (I4.16P as 

H? = (5Vj0[J + ^ o a | s=o -H/ o a + (Ug-Ug| s=o )-V||/)uS s 

-(US • V,|)Ug - (V||0° • V„)US - (Ug • V||)V,[08- (B.i) 
Here 0g = O | 2=0 - With the help of this can be further simplified to 



it; | s vf^ + v r J\j a (s')d S '\v a 0s 



-( u o • V„)Ug - (V||0g • V„)Ug - (U a • V||)V||0°. (B.2) 

Equation flM]), and the fact that u r = Vy^o imply that e 2 • curlUg = 0. Therefore, (|B~2l) can 
be rewritten as 

H? = (sV^° + V,| • I US(a')d a ' j Ug s - V,| (H£ + Ug • V„$) . (B.3) 

It follows from (I4TT8I) that 

00 00 00 s" 00 

U?| s=Q =-i y y HJ(s // )ds ,/ ds , = -^ y H?(s") y ds'cis" = -i y s"B.t(s")ds". (B.4) 

s' 
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On substituting ()B.3p into ()B.4p . we obtain 



OO OO / v. 

^ ' 



~ J sU a 0s (s) ( V|, • | Ug( S ')^'J da. (B.5) 

It follows from (jET9|) that 

00 OO 

1 sV a ds = -i^nfr a n e mT j se'^ds = -i^ner a n \ e inT = -uJ^K 



Hence, 



\ Vj0[j J s Ug (fa = -2 (V|| • <r*) a a = 2 (Vy • a a ) a%. (B.6) 

o 

Similar calculations yield 



J ^Y^sds = ^\n\\a a n \ 2 , V^- J\J%sds = -2v = 0, (B.7) 



J s U a 0s (s) ( V,| • j Ug(s') (fa' j (fa = v \n\ Re [< (V,, ■*£)]-- (V|| ■ <r a ) (B.8) 
o V o / n 

Substituting (lRBV(lR9l) into (lR5l) . we find that 

u? | s=0 = | (v,| • ^ + \ MKI a + E M Re K (V|| • K)} • (B.9) 

n n 

Finally, substitution of (|B.9|) into (|4,24p results in boundary condition (|5,2ip . 
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